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SUPERCRITICAL FLOW ABOUT A 


THICK CIRCULAR- ARC AIRFOIL 
John B. McDevitt 
Ames Research Center 

SUMMARY 


The supercritical flow about a biconvex circular-arc airfoil is being 
thoroughly documented at Ames Research Center in order to provide experimental 
test cases suitable for guiding and evaluating current and future computer 
codes. This paper is a continuation of an earlier study and extends the orig- 
inal test program to include effects of angle of attack, effects of leading 
and trailing-edge splitter plates, additional unsteady pressure fluctuation 
(buffeting) measurements and flow-field shadowgraphs, and application of an 
oil-film technique to display separated-wake streamlines. Detailed compari- 
sons of computed and measured pressure distributions for steady and unsteady „ 
flows, using a recent computer code representative of current methodology, are 
included. It was found that the numerical solutions are often fundamentally 
incorrect in that only !, strong n (shock— polar terminology) shocks are captured, 
whereas experimentally, both strong and weak shock waves appear. 


INTRODUCTION 


Rapid advances in computer technology and numerical methods lend credence 
to the belief that the computer might, by complementing the wind tunnel, 
achieve significant time and cost reductions in the design of new aircraft 
(ref. 1). Appropriate methods, however, must be devised for modeling the tur- 
bulence of compressible flows. But progress in this area is extremely diffi- 
cult and it is important that parallel experimental studies be undertaken to 
provide necessary guidance for the development and verification of newly pro- 
posed turbulence models and computer codes. 

One phenomenon of particular interest and importance in fluid dynamics is 
the separation of boundary layers by shock waves. An experimental and numeri- 
cal study of this phenomenon, as related to airfoils at supercritical speeds, 
is described in references 2 through A; the present study is a continuation of 
this earlier work. The basic approach used involved the construction of a new 
facility (the Ames "High Reynolds Number Channel"), with emphasis on achieving 
as wide a range in test Reynolds numbers as possible. This approach included 
a design concept of interchangeable test sections that become an integral part 
of each experiment. For the first experiment, a two-dimensional, 18%-thick 
biconvex circular-arc airfoil was chosen and the test-channel walls were con- 
toured to minimize wall- interference effects (see ref. 2 for details). The 
test program successfully achieved essentially two-dimensional flows over a 
wide range in test Reynolds numbers and three basic flows of interest for 


airfoils at supercritical speeds were documented: (1) steady supercritical 

flow with boundary-layer separation at the trailing-edge region, (2) steady 
flow with shock-induced separation, and (3) strong buffeting involving peri- 
odic flowfield oscillations. 

These early results, which appear to be ideally suited for guiding and 
evaluating the development of turbulence models and computer codes, prompted 
the additional studies of references 5 through 8, and gave rise to the new 
results presented here. An important modification to Deiwert's original air- 
foil code was made by Levy (ref. 6) who extended the code to include test- 
section solid-wall boundaries, thus facilitating comparisons of computation 
with experiment. Measurements of skin friction and flow-field velocities are 
presented in references 7 and 8. 

The present paper extends the experimental investigations of reference 2 
to include effects of angle of attack, effects of leading and trailing-edge 
splitter plates, additional unsteady pressure-fluctuation (buffeting) measure- 
ments and flow-field shadowgraphs, and application of an oil-film technique to 
display separated-wake streamlines. Also included here are detailed compari- 
sons of numerical solutions with experiment which will show that two fundamen- 
tal problems with numerical methods must be resolved before critical evalua- 
tions of turbulence models can proceed. 

EXPERIMENT 


The basic considerations underlying the choice of airfoil and nominal test 
Mach number are illustrated in figure 1 and the test section arrangement is 
shown in figure 2. In reference 2 a detailed discussion of the basic consider- 
ations involved in defining this experiment is presented. The airfoil model 
described in reference 2 was also used for the present tests except that addi- 
tional fast-response pressure transducers were installed for better resolution 
of unsteady-flow phenomena. The upper and lower channel wall contours, listed 
in table 1 of reference 2, were also retained for the present tests. These 
contours approximate the airfoil streamlines in free air at a nominal Mach num- 
number of 0.775. A brief summary of the steady-flow characteristics of the 
airfoil at slightly supercritical speeds will be presented before introducing 
new test results and analyses. 


Steady Flows, a = 0° 

The shock-buffeting domains for the airfoil at a = 0 , in both acceler- 
ating and decelerating flows, are shown in figure 3. At speeds below the^ 
Unsteady domain, the boundary layer separates near the trailing edge, typical 
airfoil pressure distributions are shown in figure 4. A spark-illuminated 
shadowgraph (fig. 5), taken at a moderately high Reynolds number (boundary- 
layer transition well forward of the midchord), shows multiple-compression 
wavelets that, although crowded together, do not separate the boundary layer. 

At speeds above the unsteady domains, shock- induced separation of the 
boundary layer prevails (figs. 6 and 7). The shadowgraph was taken for a 
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well-established turbulent boundary layer; a single, strong shock wave is 
clearly defined which interacts with the boundary layer to cause a large sepa- 
ration bubble with reattachment downstream of the airfoil. 


Wake Streamlines, Shock- Induced Separation 

An oil-film technique was used with thin vertical fences to display 
streamline patterns downstream of the airfoil trailing edge. Sample test 
results are shown in figure 8. The oil (a mixture of vacuum pump oil, tita- 
nium oxide, and a small amount of oleic acid) appears white against the black 
finish on the fences. The most effective display of streamlines occurred when 
the oil was applied in thin vertical lines, which served as reservoirs to con- 
tinuously feed random filaments during the test run. It should be emphasized 
that this particular application is not meaningful ahead of the airfoil trail- 
ing edge because of flow distortion in the comer between fence and airfoil, 
but downstream of the trailing edge the results are useful (see fig. 9). From 
these observations, and by making use of the LDV measurements of reference 8, 
the mean-wake— streamline pattern for "steady" flow with shock— induced separa- 
tion is approximately as that shown in figure 10. The wake is highly turbu- 
lent but the mean vertical height is approximately the same as the airfoil 
thickness and the wake closes downstream near x/c - 1.26. 

Unsteady Flow, a = 0° 

The Reynolds-number , Mach-number domains where unsteady, periodic flows 
occur for the subject airfoil was shown previously in figure 3. This unsteady 
phenomenon was discussed in considerable detail in reference 2. Additional 
information will be provided here, primarily for a chord Reynolds number of 
11 x 10 , so as to document more fully a specific unsteady-flow test case 
suitable for evaluating new airfoil codes. 

Typical instantaneous airfoil surface-pressure measurements, showing 
deviations from mean values, are presented in figures 11 through 13. The 
speed-control wedge shown in figure 2 was used to change the test Mach number 
slowly, at the approximate rate dM^/dt = 0.001. (This rate was about the 
largest that could be used without introducing noticeable errors in mean pres- 
sure measurements due to system lag.) Note that the onset of buffeting is 
sudden, but once the unsteady flow has started, the motion persists at con- 
siderably lower Mach numbers when the free-stream Mach number is decreased. 
With increasing Mach number, the phenomenon weakens and eventually steady flow 
again prevails. 

A close look at the instantaneous pressure measurements and at high-speed 
shadowgraph movies of the flow fields indicates that two distinct types of 
periodic flows occur for this airfoil at zero angle of attack. We will use 
the labels given by Tij deman (refs. 9 and 10) for various shock motions* the 
shock oscillates on the airfoil (type A); the shock disappears during part of 
the cycle (type B) ; and the shock propagates upstream and leaves the airfoil 
(type C). In the present experiment, type C motion prevails throughout most 
of the domains shown in figure 3, the exception being near the right-hand 
(i.e., high speed) boundaries where the buffeting is weak and is approaching 
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the steady state with shock- induced separation of the boundary layer. Also, 
in this case, type C motion involves alternate trailing-edge and shock- induced 
separations of the. boundary layer with the airfoil upper and lower flow fields 
180° out of phase. 

Typical shadowgraphs are presented in figure 14 ; this type of unsteady 
flow is characterized by the cyclic formation, first of a multiple— shock 
structure near the airfoil trailing edge, then emergence of a single strong 
wave, which progresses upstream and disappears forward of the airfoil. A 
dimensionless time parameter, "t, is used to indicate the fractional cyclic 
time where the start of each new cycle is coincident with the passage of the 
shock over the midchord of the airfoil. The chordwise location of the primary 
shock-wave interaction with the boundary layer for the case depicted in 
figure 14 is presented in the left-hand side of figure 15. 

The second type of unsteady flow, type A, occurs near the high-speed edge 
of the unsteady flow domain and involves only shock- induced boundary-layer 
separation. The shock wave does not disappear and reappear in cyclic fashion 
as for the first type but merely moves back and forth on the airfoil as illus- 
trated by the example presented in the right-hand side of figure 15. 

As mentioned in reference 2, the frequency of the unsteady flow was 
approximately 188 i 3 Hz throughout the test program. The slight frequency 
variation with Mach number for Re c m ** 11 x 10 6 is shown in figure 16. At 
this Reynolds n umb er the unsteady phenomenon may be summarized as follows: at 

Mach numbers below M ro = 0.733 the flow is steady with trailing-edge separa- 
tion; and above M m =°°0.778 the flow is again steady but with shock- induced 
boundary-layer separation fixed near x/c * 0.65. Between M ro = 0.733 and 
0.776 unsteady flow of type C occurs with the buffeting intensity diminishing 
with increasing Mach number. Near M ro = 0.776, the flow is erratic and the 
transition to type A begins, but the chordwise excursions of type-A motion 
di min ish rapidly with increasing Mach number and soon end at about M^ = 0.778 

We will now consider in detail the instantaneous pressure fluctuations 
for a test case involving strong, type-C buffeting at Re ro = 11 x 10 and 
a free— stream Mach number of 0.76. The cyclic variations in shock location 
and shock motion are presented in figure 17; measurements (from shadowgraph 
movies) at the shear-layer edge of shock inclination angle, 0, and of shock- 
induced flow deflection angle, 6, are presented in figure 18. Both instanta- 
neous and mean chordwise pressure variations are shown in figure 19. Pressure 
fluctuation measurements (referenced to mean values) were made at stations 
x/c = 0.25, 0.50, and 0.78 and interpolation of the data to obtain values at 
intervening stations was made with the help of information gleaned from shad- 
owgraph movies (see figs. 17 and 18). The presentation in this figure, inci- 
dentally, is arranged so that upper- and lower-surface measurements can be 
compared directly, since equivalent flow fields differ by At — 0.5. An 
alternative presentation for fixed chordwise stations is made in figure 20. 

It is interesting to note (fig. 19) that, for stations ahead of the midchord 
(airfoil crest), the maximum instantaneous pressure occurs at passage of the 
shock wave, followed by a decrease in pressure until the next cyclic shock 
arrives. However, aft of the midchord (except very near the trailing edge) 
flow compression continues after passage of the wave and the maximum pressure 
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recovery is not reached until the end of the cycle. The local flow Mach 
number, Mi, at the base of the shock wave, and the relative Mach number (shock 
fixed coordinates), 


Mi - Mi 


^shock 

a l 


are presented in figure 21. 

The information presented in figure 14 and in figures 17 through 21 con- 
stitutes an unsteady-flow test case for evaluating numerical flow simulation 
codes. A comparison of these experimental results with a current airfoil com- 
putational method will be made later in this report. 


Tests at Angle of Attack 

A limited,, and largely exploratory, study of the effects of angle of 
attack was made by rotation of the airfoil about the midchord. The channel- 
wall contours (fig. 2) were not changed and the angles of attack were limited 
to ±4°. The tests were conducted by slowly varying (dM/dt ~ 0.001) the Mach 
number from 0.750 to 0.785 and then back again to 0.750 (by using the speed- 
control wedge shown in fig. 2). The test Reynolds number was held constant at 
11 x 10 6 . From continuous recordings of instantaneous pressure measurements 
and from high-speed shadowgraph movies, the domains for unsteady flow were 
obtained and are shown in figure 22. Angle of attack serves to accelerate the 
leeward flow, enhancing the onset of shock- induced separation, but the wind- 
ward flow is retarded making it more difficult to experience shock- induced 
separation. The result is a narrowing of the' gap between left-hand boundaries 
of the unsteady flow domains (less hysteresis) as shown in figure 22. 

Typical pressure fluctuations at midchord, upper surface are presented in 
figures 23 and 24 (separate test runs were made for positive and negative 
angles of attack). The magnitude of the pressure fluctuations on the windward 
side (negative a) decrease with angle of attack. Because an angle of attack 
of 4° is relatively small for such a thick airfoil it is probable that a con- 
siderably larger angle of attack would be required for communication between 
upper and lower flow fields to weaken to where shock oscillations of type C 
could no longer be sustained. The measured frequencies of the unsteady flow 
are shown in figure 25. 

Measurements of mean pressure distributions are presented in figure 26 
for steady flow with trailing-edge separation; in figures 27 and 28 for 
unsteady flow; and in figure 29 for steady flow with shock- induced separation. 
The peak Mach number at the onset of unsteady flow remained approximately 
equal to 1.25 as the angle of attack increased but the ability to sustain 
oscillations at lower Mach numbers, once the phenomenon started, diminished 
with angle of attack. 

A composite shadowgraph (photos taken at a = 3° and -3°) is presented in 
figure 30, together with the variation of shock-wave position at the shear 
layer edge with angle of attack. 
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Splitter Plates 


The cyclic unsteady flow (type C) for the subject airfoil involves alter- 
nate trailing-edge and shock- induced separations of the boundary layer, with 
upper and lower flow fields 180° out of phase; it depends on communication 
between upper and lower flow fields. An experimental study of the importance 
of flow-field communication was made by attaching thin "splitter" plates 
(0.02 chords thick) mounted in the plane of the wing, extending upstream 
and/or downstream of the airfoil. 

The first tests were made with plates of length c/2 and typical measure- 
ments of pressure fluctuations (p. „ „ - p ) at midchord are pre- 

sented m figures 31 through 33. Mounting the plate ahead of the wing leading 
edge had little effect on the unsteady phenomenon, including flow-oscillation 
frequency (f * 190 Hz), except that the onset of unsteady flow was delayed to 
a slightly higher free- stream Mach number. As an example, for Re c = 8 x 10 6 
and dM^/dt = 0.001, the onset occurs at about M^ = 0.760 for the’wing alone 
but is delayed to about M ro = 0.765 when the splitter plate is mounted ahead 
of the wing. However, once the periodic flow oscillation starts, the motion 
persists at considerably lower Mach numbers, as in the wing-alone case. 

Attaching the plate downstream (fig. 32) virtually eliminated the 
unsteady phenomenon; the only noticeable periodic pressure fluctuations were 
detected by the most rearward sensor and this type-A motion occurred only in a 
very narrow Mach number range - from M^ = 0.771 to 0.773. When both forward 
snd rearward plates were attached, continuous monitoring of the pressure sen- 
sors in the test range from M^ = 0.745 to 0.785 failed to uncover an 
unsteady phenomenon (sample data records are shown in fig. 33). 

Since a relatively long plate in the forward location was found to have 
only a small effect on the unsteady flow, it was felt that a continuation of 
the study should be confined to downstream plate positions. Typical results 
with short trailing-edge plates are presented in figures 34 and 35 and the 
unsteady-flow domains for Re c = 8 x 10 6 are shown in figure 36. Type-C 
motion, which is easily detected by the pressure signature at x/c =0.5 (com- 
pare figs. 35 and 36), occurs only for very short plates. The effects of 
plate length on pressure-fluctuation amplitude and frequency of motion are 
summarized in figure 37. Also included in this figure is an indication of 
the downstream extent of the separation bubble for steady shock- induced sepa- 
ration at a slightly higher Mach number (from fig. 10). The wake closes about 
a quarter chord downstream of the airfoil trailing edge and a splitter plate 
of this length is very effective in suppressing communication across the wake, 
although a high-frequency shock instability of type A persists at low 
amplitude. 

The effects of various splitter plates on mean pressure measurements are 
shown in figure 38. Detailed measurements at the most rearward station (C 
at x/c = 0.975) for both increasing and decreasing speeds are presented in^ 
figure 39. Note that the large hysteresis effect in the wing-alone aft pres- 
sure measurements is greatly reduced by even the smallest splitter plate 
tested (L t e = c/8). 
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Typical spark- illuminated shadowgraphs, taken with a trailing-edge 
splitter plate of length c/2 attached, are presented in figure 40. With the 
plate, unsteady flow is avoided, except in the very narrow Mach-number range 
of h! = 0.771 to 0.773 (see fig. 36). At the lower supercritical Mach num- 
bers, multiple wavelets appear and the boundary layer is believed to separate 
only near the trailing edge. As the Mach number increases, the waves coalesce, 
first near the body, and eventually an oblique shock emerges and separates the 
flow. 

COMPARISONS OF NUMERICAL SOLUTIONS WITH EXPERIMENT 


Three supercritical flow fields will be examined in detail: (1) the 

supersonic zone terminated by diffuse pressure recovery with attached boundary 
layer; (2) shock- induced separation of the boundary layer; and (3) unsteady 
asymmetric flow involving periodic switching between trailing-edge and shock- 
induced separations. Experimental measurements will be compared with numeri- 
cal solutions from an airfoil computational method developed by Deiwert 
(ref. 3), which includes MacCormack's modifications (ref. 11) to improve effi- 
ciency, and Levy*s modification (ref. 6) to include test-section solid-wall 
boundaries. Basically the method solves the time-dependent, two-dimensional, 
Reynolds-averaged Navier-Stokes equations with simple algebraic eddy viscosity 
models (see refs. 2 through 6 and 11 for further details). 


Trailing-Edge Separation 

In the Mach-number range from critical (first appearance of sonic flow on 
the airfoil) to the onset of buffeting, the airfoil supersonic zones, with 
adjacent turbulent boundary layers, terminate by diffuse pressure recovery 
and the boundary layer separation is delayed to near the trailing edge. 

When a trailing-edge splitter plate is used, the unsteady flow is suppressed 
and the shadowgraphs presented in figure 40, all for steady flow, show the 
transition from multiple shocks to a single oblique shock as the Mach number 
increases. 

A comparison between experimental and computed pressure distributions for 
the case where multiple wavelets appear, but do not separate the boundary 
layer, is presented in figure 41. The comparison is made for the wing with 
trailing-edge plate included because Deiwert 1 s code indicated unsteady flow 
for the wing alone at this Mach number (experimentally, either steady or 
unsteady flow can occur; see fig. 3). The measurements indicate a steep pres- 
sure recovery between about x/c = 0.6 and 0.68 where the wavelets are most 
concentrated and the pressure recovery is matched remarkably well by computa- 
tion. Although the computational grid is closely spaced vertically near the 
airfoil surface to resolve the boundary layer, the horizontal spacing is rela- 
tively large (Ax/c - 0.02 near the "captured shock" and larger elsewhere) - 
this tends to smear the pressure recovery but there was good agreement between 
experiment and computation except very near the trailing edge where the 
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boundary layer had separated. The poor agreement near the trailing edge 
suggests inadequate turbulence modeling. 

The complex flow phenomenon indicated by the shadowgraph, rules out "shock- 
fitting" techniques. The "shock-capturing" technique (included i n the Deiwert 
code) probably is sensitive to grid geometry in many cases and a definitive 
study of this aspect would be particularly informative using the complex flows 
displayed in the photographs of figure 40 as examples. 

The computed Mach contours in the flow field are shown in figure 42. The 
contours stand more or less normal to the surface near the downstream edge of 
the supersonic zone (x/c * 0.65) but the rapid pressure recovery is not repre- 
sentative of a single normal shock. Specifically, the maximum flow Mach 
number, both the computational number and the experimental number, is about 
1.26 and the pressure jumps from about Cp = —1.1 to — 0.4j the pressure jump 
for normal shock recovery, however, would be from -1.1 to -0.2. 


Shock— Induced Separation 

A comparison of measured and computed pressure distributions for steady 
flow at M^ = 0.785 with shock- induced separation is presented in figure 43. 

The numerical method accurately predicts the pressure distribution over the 
forward portion of the airfoil but is grossly inaccurate in predicting shock 
location and strength. The fundamental differences between the numerical and 
experimental shock waves will be discussed next. 

Computed Mach contours at M ra = 0.783 are shown in figure 44(c) (repro- 
duced from ref. 6). In figure 44(c), "chords traveled" is a measure of elapsed 
computational time after convergence of solutions and corresponds to chord 
lengths traveled by the mean flow during one cycle of the oscillatory flow - 
figure 44(b) - to be discussed later. The numerical solutions at = 0.783 
behaved erratically, but a periodic unsteady flow was not indicated, and in 
the following analysis mean values of computed solutions at this Mach number 
will be used. (The reason for this erratic behavior in the solutions at this 
Mach number will be discussed subsequently.) In this flow example an oblique, 
clearly defined shock wave appears and the pressure recovery near the foot of 
the wave is compatible, both experimentally and computationally, with theoreti- 
cal oblique shock relationships, as shown in figure 45. In figure 45, com- 
puted and measured values of shock wave angle, pressure jump, and local flow 
Mach number ahead of the shock have been superimposed on the theoretical 
inviscid relationships. The experimental measurements locate -a point in the 
"weak" shock-polar domain and the flow downstream of the wave is supersonic. 

The computed values were taken to be the mean values of the somewhat erratic 
solutions shown in figure 44(c) and define a location in the "strong" shock- 
polar domain for which the downstream flow is subsonic. 

The more conventional ways to display shock-polar relationships involve 
the flow deflection angle, as shown in figures 46 and 47. The points nearest 
the airfoil, "A" (experiment) and "E" (computation), were transcribed from the 
previous figure and imply approximately the same flow deflection angle in both 
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cases (6 * 7°). Sketches of the computed and experimental flow fields, drawn 
to scale, are included in figure 46, In the computation, the shock wave is 
always of the "strong type," supersonic to subsonic, and numerical irregulari- 
ties in the subsonic downstream region can directly affect the shock wave. 

This is the probable explanation for the quasi-steady behavior mentioned pre- 
viously of the computations shown in figure 44(c). Experimentally, the shock 
is "strong" only in the outer region, (C) to (D) , while inboard the shock is 
"weak" with the so-called "supersonic tongue" attached from (A) to (B) . This 
"tongue" acts as a buffer region between the turbulent wake and shock wave 
which explains why the shock position appears to be very stable in the high- 
speed shadowgraph movies at this test condition. 

Shock-polar variations between downstream Mach number and flow-deflection 
angle, for various values of upstream Mach number, are shown in figure 47 (the 
flow-field points (A) through (E) correspond to those shown in the previous 
figure). Measurements by a laser Doppler velocimeter technique (ref. 8) of 
local Mach number at the shear-layer edge are presented in figure 48. These 
measurements indicate that the "supersonic tongue" extends nearly to the air- 
foil trailing edge. 

The poor agreement between experiment and computation shown in figure 43 
is a consequence of the numerical method choosing a "strong" (shock-polar ter- 
minology) shock wave whereas the experimental wave is primarily "weak." It is 
well known from shock-polar relationships that two different shock waves and 
flow patterns are theoretically possible for given 6 and Mi . However, it is 
believed that only the weaker shock wave (larger M 2 and smaller 0) will 
occur experimentally (see discussions of flows past wedges in ref. 12, p. 12, 
and ref. 13, p. 317). It is doubtful that computational grid size or orienta- 
tion could account for the fact that the numerical method captures the shock 
on the wrong side of the shock polar. 


Unsteady Flow 

Comparisons of measured and computed pressure distributions at a Mach 
number where strong buffeting occurs are presented in figure 49. (See figs. 17 
through 21 for additional experimental results at this Mach number.) Deiwert’s 
numerical method (ref. 3) has been shown by Levy (ref. 6) to be capable of 
qualitatively predicting the essential features of this unsteady flow but 
significant differences between experiment and computation are evident, par- 
ticularly with regard to the aft surface pressure recovery. The frequency of 
pressure oscillations was also noticeably different — in the experiment the 
reduced frequency was f = 0.49 (f ~ 190) but numerically the reduced frequency 
was about 20% less (f ~ 150) . In both the experiment and in the numerical 
solutions a clearly defined shock wave emerges after an elapsed cyclic frac- 
tional time of about "t = 0.25. Detailed comparisons of shock positions at 
the shear-layer edge and of corresponding effective upstream Mach numbers and 
shock pressure jumps are presented in figures 50 and 51. In both the experi- 
ment and computations the cyclic-shock formation occurs near x/c = 0.85, but 
in the experiment the shock wave immediately after formation moves upstream 
rapidly, whereas the computed shock wave remains stationary near x/c - 0.85 
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for about 20% of the cycle. This difference in shock-position behavior 
probably is related to the fact Cto be shown subsequently) that the shock wave 
near the shear-layer edge experimentally is "weak" Cshock-polar terminology), 
with trailing flow supersonic, whereas the computed shock is "strong," with 
trailing flow subsonic. Consequently, the stream momentum external to the 
separated region is much higher experimentally; this supports a stronger 
cyclic- induced upwash/downwash phenomenon and the experimental flow oscilla- 
tion does not hesitate as in the numerical solutions. 

Experimental measurements at or near the airfoil surface of (p 2 - Pl)/<ll> 
Mi, and Q at the shear-layer edge are superimposed on shock polar diagrams in 
figure 52 and these three measurements are, within the accuracy of the tests 
and for the cyclic times considered, mutually compatible with pressure recov- 
ery through oblique shock waves. The results are summarized in figure 53 and 
compared with numerical solutions. In the computations, only "strong shocks" 
appear, but in the experiment the first cyclic appearance of a clearly defined 
wave is in the "weak-shock" domain of the shock polar. With increasing cyclic 
time, the strong-shock domain is entered as the wave approaches the airfoil 
crest and eventually disappears ahead of the airfoil as a weakening sound 
wave. The failure of the numerical method to capture anything but strong 
shocks is the probable explanation for the inaccurate pressure recoveries 
shown in figure 52. 

The unsteady flow analyzed here has not introduced new computational 
problems. The failure of the numerical code in the unsteady case is related 
to the failure discussed previously for steady flow with shock- induced separa- 
tion in that only "strong" shocks are found in the computations. If a suit- 
able computing code and adequate turbulence modeling are found for the static 
case, the dynamic case should, in principle, follow, as time-accurate solu- 
tions are now routinely obtained with advanced computers. 


CONCLUDING REMARKS 


The present study of the supercritical flow about a thick circular-arc 
airfoil is a continuation of the investigation reported on in reference 2. 

The primary purpose of this undertaking is to provide well-documented experi- 
mental test cases suitable for guiding the development of and evaluating air- 
foil numerical flow- simulation codes. This paper extends the original test 
program to include effects of angle of attack, effects of leading and trailing- 
edge splitter plates, additional unsteady-flow pressure measurements and flow- 
field shadowgraphs, and application of an oil-film technique to display 
separated-wake streamlines. 

Included in this paper are detailed comparisons of computed and measured 
airfoil surface pressures for three specific test cases, all involving super- 
critical flows. The numerical method used in the comparisons represents cur- 
rent methodology and solves the Reynolds-averaged, two-dimensional, compressi- 
ble Navier-Stokes equations. The first test case studied involved steady flow 
with diffuse pressure recovery terminating the supersonic zone and with 
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boundary layer separation delayed to near the trailing edge. With the 
exception of the trailing-edge region, the agreement between computed and 
measured pressures was good. The poor agreement for the separated trailing- 
edge region indicates inadequate turbulence modeling. 

The next test case analyzed involved shock- induced separation of the 
boundary layer. Here, the numerical simulation was found to be fundamentally 
incorrect in that the method captured a "strong shock" (shock-polar terminol- 
ogy) whereas the experiment showed a "weak shock," with trailing "supersonic 
tongue" near the body, and the shock wave is "strong" only near the outer 
extent of the supersonic zone. Consequently, the computed aft surface pres- 
sure recovery was much larger than that measured. 

The third case involved strong buffeting with separation alternately 
switching from trailing-edge separation to shock- induced separation. Again, 
it was found that the numerical method captures only strong shocks and that 
the experiment displays strong and weak waves in the cyclic-flow oscillation. 

The detailed comparisons of computed and measured results point out the 
need for two fundamental studies in addition to the general problem of finding 
proper turbulence models: first, the influence of computational mesh geometry 

on the numerical solutions must be determined. Second, the preference of 
numerical methods for strong shocks, in contrast to the appearance experimen- 
tally of both weak and strong shocks, must be corrected. Unsteady flow phe- 
nomena, such as the buffeting included in the present study, do not introduce 
additional computational problems. If the numerical method and the turbulence 
modeling are adequate in the static case, then time-dependent solutions would 
easily follow. 
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- Pressure distributions at speeds below the buffeting region 

(from ref. 2). 


16 






0 .2 .4 .6 .8 1 

x/c 


Figure 6.- Pressure distributions for steady flows at speeds above the buffet- 
ing region (from ref. 2). 
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Figure 7.- Flow-field shadowgraph, = 0.785, Re c ^ = llxlO 6 . 
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Figure 8.- Oil-flow patterns on vertical fences 
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Figure 10.- Wake streamlines, shock- induced separation; = 0.785, 

Re Cj0O = UxlO 6 . 
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(b) Re c ^ = 5xl0 6 to 8xl0 6 . 


Figure 13.- Concluded. 












Figure 14.- Unsteady flow-field shadowgraphs; = 0.76, Re 








' : ' V '-V. V< - - s*-, -- ’ 


%r <j.Ay /' ,C' ' 

Afciffi 






M 


lliil 




Continued 









mm. 




Continued 


















Concluded 







Figure 15.- Shock motions in unsteady flow, Re c m = 11x10®. 
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Figure 17.- Cyclic shock motion; M ra = 0.76, Re Cj0o - 
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Figure 18.- Cyclic shock wave and flow deflection angles; = 0.76, 

Re_ = llxlO 6 . 
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Figure 19.- Mean and instantaneous airfoil chordwise pressure distributions in 
unsteady flow; = 0.76, Re m = 11*10 6 , f = 189 Hz (T TO ~ 270 K) . 















Figure 20.- Time histories of airfoil pressures; Me = 0.76, Re c m - llxlO 6 , 

f = 189 Hz. 
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Figure 24.- Effect of Reynolds number on pressure fluctuations at a — 4°, 

M*, = 0.76. 


41 









Figure 26.- Airfoil pressure distributi 
with trailing-edge separation; 


at angle of attack for steady flow 
= 0.745, Re„ = 11x10 s . 
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Figure 30.- Composite shadowgraph of flow field at a = 3 
shock location with angle of attack; If “ 0.785, Re 


and variation of 
= 11x10® . 
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Figure 35.- Pressure fluctuations with trailing-edge plate of length 
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(a) Wing alone. 

Figure 38.- Mean pressure distributions with various splitter plates 






(b) Leading-edge plat 
Figure 38.- Continued 
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Figure 39.-* Variation of airfoil pressui 
Mach number for variot 
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Flow-field shadowgraphs of airfoil with trailing-edge plate of 
length c/2; Re = 11x10^. 
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Continued 
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Figure 40.- Continued, 
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Figure 43.- Comparison of computation and experiment; = 0.78, Re c ^ = 10 7 . 
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STEADY FLOW, TRAI LING-EDGE SEPARATION, = 0.720 
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Figure 45.- Oblique shock solutions, (p£ - pjJ/q^ vs 0 for given , and 
overlays of computed and experimental values for = 0.785 and 
Re = llxlO 6 . 
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Figure 46.- Oblique shock solutions, 9 vs <S for given M, , and overlays of 
computed and measured values for = 0.785 and Re* m = llxlO 6 . 
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Figure 53.- Oblique shock relationships and comparison of computed and experi- 
mental measurements at the shear-layer edge in unsteady flow; = 0.76, 

Re = llxlO 6 . 

C y OO 


79 



2. Government Accession No. 


3. Recipient's Catalog No. 


1. Report No. 


4. 


NASA TM- 78549 

Title and Subtitle 


5. Report Date 


SUPERCRITICAL FLOW ABOUT A THICK 
CIRCULAR-ARC AIRFOIL 


6. Performing Organization Code 


7. Author(s) 


8. Performing Organization Report No. 


John B. McDevitt 


9. Performing Organization Name and Address 

NASA Ames Research Center 
Moffett Field, Calif. 94035 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


A- 769 3 

10. Work Unit No. 

505-06-11 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 


16. Abstract 


The supercritical flow about a biconvex circular-arc airfoil is being 
thoroughly documented at Ames Research Center in order to provide experi- 
mental test cases suitable for guiding and evaluating current and future 
computer codes. This paper is a continuation of an earlier study and 
extends the original test program to include effects of angle of attack, 
effects of leading and trailing-edge splitter plates, additional unsteady 
pressure fluctuation (buffeting) measurements and flow-field shadowgraphs, 
and application of an oil-film technique to display separated-wake stream- 
. Detailed comparisons of computed and measured pressure distribu- 
tions for steady and unsteady flows, using a recent computer code 
representative of current methodology, are included. It was found that 
the numerical solutions are often fundamentally incorrect in that only 
strong (shock-polar terminology) shocks are captured, whereas experimen- 
tally, both strong and weak shock waves appear. 


17. Key Words (Suggested by Author(s)) 

Fluid mechanics Supercritical flow 

Numerical methods Separated flow 
Airfoils Unsteady flow 

18. Distribution Statement 

Unlimited 

STAR Category - 02 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

83 

22. Price* 

$5.00 i 


*For sale by the National Technical Information Service, Springfield, Virginia 221 61 















NASA Technical Library 



01473 6426 



